Simulations on the Heavy Hadron Transport at RHIC 
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Based on the hadron transport frames, detailed simulations are carried out to discuss (j>, ^ pro- 
ductions and the significant enhancements in the very low pr region for some of the soft spectra in 
RHIC. Elastic interactions are introduced in the simulations. The elastic cross sections vary from 
different hadrons and energy scales, which can qualitatively explain the different collective motions 
of various hadrons. 
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INTRODUCTION between different hadrons and pr regions. 



In the recent progress [![, we have developed a trans- 
port model to describe the hadron production in rela- 
tivistic heavy ion collisions. In the model, a decoupling 
hypersurface of the fluid which consists of a group of 
splitting QGP droplets Q is used to emit hadrons. The 
model satisfies most of light hadrons' spectra, but dis- 
agrees with the distributions of </> [H 0] , O || and some 
heavy resonances 0,0). As another remarkable phe- 
nomenon, there seems to be significant enhancements at 
very low transverse momentum in some of the spectra 
for A 8], S*(1385) and <fi, to change the slopes. On the 
contrary, the model predicts a decrease when px drops. 

The idea of the hypersurface hadronization locks the 
emitted hadrons in the local rest frames of the hypersur- 
face. This is right when the fluid is dense enough in the 
early stage of the evolution. In the later stage, when 
the inelastic free path becomes as long as 4 /to, the 
emission fixed in the hypersurface might be less strict. 
Considering the spectra of light hadrons are well satis- 
fied, as other hydrodynamic methods with decoupling 
hypersurface did, there must be some reasons to push 
the inside light hadrons to participate in the collective 
motions. One of the possible candidate is the elastic col- 
lision j|. Although elastic collisions do not change the 
hadron abundances themselves, the effect of preventing 
hadrons from entering the more dense regions and being 
absorbed, will enhance their abundances or momenta. 

Thus, detailed simulations on hadron transport are ap- 
plied to some of the heavy hadrons, such as <f> and f2 in the 
central events in RHIC Au-Au 200 GeV . The advantages 
for <p and f2 are based on the condition that almost no 
hadrons could decay to them and their widths are large 
enough to neglect the affections of decays in the expand- 
ing medium. As a simplification, only selected hadrons 
for the measurements are simulated. The medium evolu- 
tion and the feedbacks to the medium are not simulated 
in microscopic frames. As the cross sections of some 
heavy or multi-strange hadrons may be much smaller 
than those of light hadrons [10], it will provide us a pos- 
sible solution to explain the various collective behaviors 



MEDIUM EVOLUTION 

The decoupling of the fluid is a gradual progress in 
principle. As the feedback of the hadron simulation to 
the fluid is unjustified, the medium is separated to an 
ideal hydrodynamic part and a free expanding part with 
evaporation by a cut e?,. The 1 + 1 hydrodynamic evo- 
lution with an invariant boost for the central events in 
RHIC Au-Au 200 GeV is determined by the method in 
ref [l|. The decoupling surface here is similar to the re- 
sults of surface evaporation [l[ except in the very early 
stage which is considered to be less important. For a 
group of decoupled droplets, there are no collective mo- 
tions and local rest frames indeed. Approximately, the 
whole evolution is defined, 



dl= {£cVR — 



e c n vR ^^-}dRdv, (1) 



where the former term in the integral stands for the free 
flights of the droplets and the latter term is the evapora- 
tion. Then, 
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is the velocity and radius distri- 



bution of the number density of the droplets. 

As we do not know the distribution of the droplets 
n V R and can not simulate all of the droplets cither, the 
expansion is simplified as, 
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where X — and N = 2 is the symmetric dimension. 
Considering that the droplet evolution by emitting and 
absorbing hadrons is not clear, the evaporation term is 
omitted in our calculations. 
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HADRON TRANSPORT 

The classical transport equation for the evolution of 
hadron phase space distribution /(£, r, p) can be written 
as 



|(-+ v - Vf = Q + a-Pf, 



(5) 



where the second term on the left hand side is the free 
streaming part, and on the right hand side the lose and 
gain terms /3(p)/ and a(p) indicate hadron absorption 
and production in the medium, and Q(p; f) stands for 
the elastic term. 

In the simulations, hadrons generated in QGP or 
decayed from heavier hadrons, will pass through the 
medium and decay at the same time. Thus, the gain 
term includes two parts, a = uqqp + adc- The decay 
contribution is easy to obtain, but the term (Xqgp is 
complicated. The reason is that the quantitative tem- 
perature dependencies for some parameters are unknown 
especially in a rapidly expanding system. Therefore, 
the production is simplified as ideal QGP with massless 
quarks and gluons plus the contributions of bound states 
which will never change via T, as well as the coupling 
constants. The simplification is based on such an idea 
that in large T regions, where the inelastic free path is 
short and the produced hadrons may travel more distance 
before they escape from the medium, the affections may 
be less important for the simulations. Thus, the approx- 
imation will be sufficient enough if the productions are 
well satisfied near the critical temperature T c . Then, 
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where I = 2 for mesons and I = 1 for baryons due to the 
quark-diquark frame [l2| . 77 is the ratio of the occupation 
volume of all the QGP droplets to the whole volume of 
the system, «o is the production rate at T c and 
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where tBound is contributions of bound states in QGP. 

The loss rate includes the decay rate and the medium 
absorption, (3 = (3 dc + f3 abs , where f} dc = 7 _1 r and 
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(7) 



with 0o w {\ + v 2 /"i)n q {u) i ne i as tic at the critical tempera- 
ture T c . As the expression belongs in the local rest frame 
of the medium, a Lorentz transformation on time step is 
required when the simulation is applied in the selected 
frame. Details on a and Po are discussed in ref [J, [13] • 
It is noted, © and (J7JI are not independent, even the 
droplets are assumed under random distribution. When 



77 < 1, the location of production and absorption is cor- 
related. If the QGP droplets are large enough, a factor 
k(R) should be inserted to the independent productions 
([5]) to remove the absorption by the droplet where the 
hadron is produced by QGP combination in principle. 
k(R) can be fitted as, 
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where 
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(9) 



As discussed in above sections, the evolutions of the 
droplets are unjustified so far. Thus we set R ~ and 
k(R) — 1 to minimize the parameter dependency. 

Hadrons produced in time, position and momentum 
bins (?) are tracked when they pass through the medium. 
They are recorded separately, as 



Ni = N escape + N r 



(10) 



where N escape is the direct and reconstructed hadrons 
escaped from the medium by cuts and N reco is the re- 
constructed hadrons which may decay in the medium, 
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(ii) 



with recorded decay products rij and 

N[ = Decayed hadrons x Branch ratio. (12) 

There is a small difference between our calculation 
and the experimental manipulations. In the experiments, 
usually one channel is measured for the reconstruction to 
provide the total abundance by dividing the correspond- 
ing branch ratio. 

Inelastic simulations are applied with elastic cross sec- 
tions omitted (Q = 0). For the consistency with [H, no 
parameter is changed. There is no new parameter ei- 
ther. The results for and f2 are shown in Fig. [T^, and 
Fig. O The simulated results look much better than 
the original estimations from hypersurface emissions [l[. 
The productions in very low pt region are enhanced and 
the collective flows are reduced. The perfect f2 spec- 
trum shows the elastic collisions for could be neglected 
as they may not participate in the collective motions. 
As the enhancement in very low pt region for <fi is still 
smaller than expected, elastic collisions are required in 
our considerations. 



ELASTIC COLLISIONS 

The elastic process can be divided to two effects ap- 
proximately. The acceleration pushed by the pressure 
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of the medium, 
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FIG. 1: Model estimations for mr distribution of (A(1020) 
comparing with the data from STAR Collaboration 0] and 
the PHENIX Collaboration (a) Comparison between hy- 
persurface hadronization with transport emit distributions |l| 
and inelastic transport simulation, (b) Demonstration of the 
push and drag effects by elastic collisions, (c) Demonstra- 
tion of drag effects at different mean elastic cross sections. 
Subscripts for the mean elastic cross sections are removed for 
clarity. 
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FIG. 2: Model estimations for pr distribution of Q + Q + 
comparing with the data from STAR Collaboration [5j. 



gradient and the deceleration dragged by the relative mo- 
tion. Different from the inelastic cross sections, which is 
nearly stable via energy scales and parton species, the 
elastic cross sections vary from 10 1 mb to 10 3 mb [l3T |. 
which makes it difficult to obtain a strict calculation with 
limited information of parton-hadron or hadron-hadron 
elastic differential cross sections, especially in the condi- 
tion that only a few selected hadrons are simulated. 

Heavy hadrons are discussed here to present a rough 
estimation on elastic process. As they are heavy, the 
changes of their momenta after each collision and the 
thermalization are supposed to be small. Thus, the push 
and drag effects can be estimated by integrating all of 
the collisions or the probabilities of them. 

The push effect can be described in the local rest frame 



dpi 

dt 



(13) 



push 



where (V)i is the mean hadron volume of elastic colli- 
sions. As detailed differential cross sections are not clear, 
it is estimated by a rigid ball approximation, 
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On the other hand, the drag effect for a flat plate mov- 
ing in normal direction is presented as, 
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Jd 3 qjfj(qj)\vj - Vi||| Apijipi,^), 



where Apjj(pj, qj) is the momentum shift after each col- 
lision. The integral could be simplified by assuming the 
medium particles to be massless and fixing some of the 
momentum symbols to a kind of average (pj). Thus, 
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where, 



E l (E l + (p J )) 

[Ei(l + v t ) + - Vi ) + fa)]' 



(14) 



(15) 



with (pj) « 500 MeV for medium particles near T c ss 
166 MeV . Cf is a correction function for the "shape" of 
the hadrons. 

It should be noted, as the elastic cross sections vary in 
a large range for different energy scales, different averages 
of elastic cross sections (cr) ip and (cr}id will be reached by 
equations (fTB")) and Q3]). The relation between them is 
not clear. As Eq p3[) affects slow hadrons more inten- 
sively and Eq (jT4j) affects fast hadrons more significantly, 
(a)ip might be slightly larger than {a) id- To make the 
value (a}id similar to (a)i P , we estimate c/ ~ 0.25 as a 
rough approximation. 

The elastic corrections are only applied in the medium 
before it decouples. Equations (fTB")) and (fT4"|) are inserted 
to the simulation separately. The results are troublesome 
when the push and drag effects are inserted together. 
They disturb each other in their invalid regions and the 
problem is not resolved so far. As shown in Fig. [JJd, 
4> production at very low pr is enhanced by the drag 
effect of preventing hadrons from moving inside. The 
enhancement 0] could be well satisfied by fitting {a) ic i — 
20 mb. In higher px regions, push effect will effectively 
accelerate the hadrons, like the hydrodynamics does, to 
participate in the collective motions until the medium 
decouples, when {<j)i P = 20 ~ 25 mb. They don't work 
well in opposite regions. The effects of different {o)id are 
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Although the equations (| 13[) and (fT4)) are not suitable 
for light hadrons, we still applied them to test the push 
effect for pions. When {a) ip ~ 75 mb, the inclusive pion 
distribution 14J can be reproduced, as shown in Fig. 
HI For direct pions without any decay contributions, 
(<?) 



250 mb. Local thermalization as a Boltzmann 
distribution was tried for strongly interacting regions for 
all the simulations and no significant difference or better 
results were found. 



SUMMARY 



FIG. 3: Model estimations for inclusive mr distribution of 
(1385) + E* T (1385) comparing with the data from STAR 
Collaboration [(j}. {a} for mesons is (|) 2,/3 of the values listed. 



shown in Fig. [Tfc. The result for implies the elastic 
cross section for O is no more than 20 mb. 

As the momentum of each surviving daughter hadron 
may be shifted by elastic collisions, it is hard in cur- 
rent simulations to determine the efficiency of the recon- 
struction and the changed momenta of the reconstructed 
hadrons, which are different from those when they decay, 
until all of the related hadrons including backgrounds 
are simulated to commit a real reconstruction. Further 
more the elastic cross sections may be much different for 
hadrons and their decay products. Thus, The spectra for 
some strong decay hadrons, can not be simulated strictly. 
The spectrum of E*(1385) is not satisfied at very low pr 
so far, as shown in Fig. [3] S*(1530) has the same prob- 
lem but looks more serious. Other possibilities may be 
the lack of information for some of heavier baryons which 
may decay to S*(1385) or S*(1530), and inelastic hadron 
re-scattering. 
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FIG. 4: Model estimations for inclusive pr distribution of 
comparing with the data from PHENIX Collaboration [l4| 
(a) for baryons is (|) 2//3 of the values listed. 



Beyond the approximation of hypersurface hadroniza- 
tion developed from the transport model we have 
applied transport simulations for heavy hadrons to repro- 
duce the significant enhancement in very low px region 
in some of transverse momentum spectra. Although the 
method is neither strict nor complete, and the values of 
mean elastic cross sections are just qualitative, the reason 
for the enhancement and therefore the different slopes in 
the spectra may be successfully explained by the push 
and drag effects. The simulation shows that the elastic 
cross sections of <p and Q are much smaller than those of 
other hadrons. The simulation can also explain the rea- 
son why the hypersurface hadronization works well for 
light hadrons by the push effect. 
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